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Abstract. This is a survey on tropical polytopes from the combinatorial point of view 
and with a focus on algorithms. Tropical convexity is interesting because it relates a 
number of combinatorial concepts including ordinary convexity, monomial ideals, sub- 
divisions of products of simplices, matroid theory, finite metric spaces, and the tropical 
Grassmannians. The relationship between these topics is explained via one running cx- 
' ample throughout the whole paper. The final section explains how the new version 2.9.4 

, of the software system polymake can be used to compute with tropical polytopes. 

(N 



00 



o 

(N 



>< 



1. Introduction 



The study of tropical convexity, also known as "max-plus convexity", has a long tra- 
dition going back at least to Vorobyev [38] and Zimmermann [41]; see also [9l [26l [8] 
O ' aiid the references there. Develin and Sturmfels contributed an inherently combinato- 
^ ■ rial view on the subject, and they established the link to tropical geometry [12]. The 
! subsequent development includes research on how tropical convexity parallels classical 
^ \ convexity [201 ID Ell [H] and linear algebra [34] as well as investigations on the relation- 
ship to commutative algebra [5l [13]. Most recently, topics in algebraic and arithmetic 
geometry came into focus [251 ISSl [221 IH]- The purpose of this paper is to list known 
^ ! algorithms in this area and to explain how the connections between the various topics 
^ I work. We emphasize the aspects of geometric combinatorics. 

Our version of the tropical semi-ring is (M, min, +), and we usually write "©" instead of 
. "min" and "0" instead of "-[-". Of course, it is just a matter of taste if one prefers "max" 
! over "min". Whenever convenient we will augment the semi-ring with the additively 
0> I neutral element oo which is absorbing with respect to 0. The tropical d-torus is the set 
§ ■ T'^"^ := M^I^ld, where 1^ := (1,1,..., 1) is the all-ones-vectors of length d. Via the 
O ■ maximum-norm on R'^ and the quotient topology the tropical torus T'^~'^ carries a natural 
topology which is homeomorphic to W^~'^ . Componentwise tropical addition and tropical 
scalar multiplication turn W'- into a semi-module, and since tropical scalar multiplication 
■ with A is the same as the ordinary addition of the vector A ■ 1^ tropical linear combinations 
^ of elements in 1"^-! are well-defined. For a set \^ C T'^"^ of generators the tropical convex 
hull is defined as 

tconv V := {XQv® fJ^Qw \ X, fi e R, v,w e V} . 

If V is finite then tconv V is a tropical polytope, which can also be called a "min-plus 
cone". There are several natural ways to represent a tropical polytope, for instance, as 
the tropical convex hull of points (as in the definition), as the union of ordinary polytopes 
|12] . or as the intersection of tropical halfspaces [2U]. We will discuss several "tropical 
convex hull algorithms", that is, algorithms which translate one representation into the 
other. This turns out to be related to algorithms in ordinary convexity as well as to 
algorithms known from combinatorial optimization. Here we focus on the key geometric 
aspects. 
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Combinatorics enters the stage through the observation [121 Theorem 1] that configura- 
tions of n points in T'^"^ or, equivalently, tropical polytopes with a fixed set of generators, 
are dual to regular subdivisions of the product of simplices A„_i x A^-i. Products of sim- 
plices in turn occur as the vertex figures of hypersimplices, and these hypersimplices are 
known to serve as adequate combinatorial models for the Grassmannians. In fact, it turns 
out that the tropical Grassmannians can be approximated in terms of decompositions of 
hypersimplices into matroid polytopes [371 136| [T7] . We will show how a configura- 
tion of n points in 1"^-! can be lifted to a matroid decomposition of the hypersimplex 
A(d,n + d). 

The structure of this paper is as follows. The first section briefiy collects some informa- 
tion about the tropical determinant and tropical hyperplanes before the next one intro- 
duces the basic combinatorial concepts of tropical convexity. This section also explains 
how ordinary convex hull algorithms can be used to compute with tropical polytopes. 
Then we discuss various versions of the tropical convex hull problem. In a section on 
matroid subdivisions we explore the role of tropical polytopes for the tropical Grassman- 
nians and finite metric spaces. Finally, we show how the new version 2.9.4 of the software 
system polymake can be used for computing with tropical polytopes. 

2. Tropical Determinants and Tropical Hyperplanes 

Evaluating ordinary determinants is the key primitive operation of many algorithms in 
ordinary convexity. Its tropical counterpart is equally fundamental. The tropical deter- 
minant of a matrix M = {rriij) G M"'^'^, also called the "min-plus permanent", is defined 
by the tropicalized Leibniz formula 

(1) tdet M := ^ rrai^^(i) m2,^(2) ■ ■ ■ md,a{d) , 

o-eSym^ 

where Sym^ is the symmetric group acting on the set [d] := {1,2, ... ,d}. Evaluating tdet 
is the same as solving the linear assignment problem, or "weighted bipartite matching 
problem", from combinatorial optimization for the weight matrix M. This can be done 
in 0{d^) time [331 Corollary 17.4b]. 

By definition det M vanishes if the minimum in the defining Equation ([T]) is attained 
at least twice. In this case M is tropically singular, and it is tropically regular otherwise. 
Checking if the tropical determinant vanishes can be translated into solving d + l assign- 
ment problems as follows: First we evaluate tdet M be solving one assignment problem. 
This way we find some optimal permutation a with tdet M = mi^o.(i)0'^2,o-(2)©' ■ ■&^d,a{d)- 
The permutation a is called a realizer of tdet M. We have to check if there are other re- 
alizers or not. To this end define d matrices Mi, M2, . . . , Md all of which differ from M in 
only one coefficient, namely the coefficient in the z-th row and the (T(i)-th column of Mj 
is increased to mj_o-(j) + 1. We subsequently compute tdet Mi, tdet M2, ... up to tdet M^. 
Clearly tdet M < tdet Mj for all i. Now M is tropically singular if there is another per- 
mutation r for which the minimum is also attained. As a and r must differ in at least 
one place it follows that such a r exists if and only if tdet M = tdet Mj for some i. We 
conclude that M is tropically regular if and only if tdet M < tdet Mj for all i. Hence 
deciding if M is tropically singular or not requires 0{d^) time. 

Remark 1. With the tropical determinant we can express that a square matrix in the 
tropical world should be considered as having full rank or not. Defining the rank of a 
general matrix in the tropical setting is much more subtle |llj . 

The tropical evaluation of a linear form a G M'^ at a vector x G T'^"^ reads (a,x)trop '■ = 
Oi Xi © 02 X2 © ■ ■ ■ © Xrf. Again this expression vanishes if the minimum is attained 
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at least twice. A vector a G M defines a tropical hyperplane 

H{a) := {x G T'^ | (a, x)trop vanishes} , 

and the point —a G T'^"^ is called the apex of H{a). Any two tropical hyperplanes just 
differ by a translation. We have the following tropical analog to the situation in ordinary 
linear algebra; see also [HI [18] . 

Proposition 2 ( [30l Lemma 5.1]). The matrix M is tropically singular if and only if 
its rows (or, equivalently, its columns) considered as points in T*^"^ are contained in a 
tropical hyperplane. 

The ordering of the real numbers allows to refine the linear algebra in M'^"^ to the 
theory of ordinary convexity. For algorithmic purposes sidedness queries of a point versus 
an affine hyperplane (spanned by d — 1 other points) are crucial, and this is computed by 
evaluating the sign of the ordinary determinant of a d x d-matrix with homogeneous coor- 
dinate vectors as its rows. This does not directly translate to the tropical situation since 
the sign of the tropical determinant does not have a geometric meaning. More severely, 
the complement of a tropical hyperplane in T'^-^ has exactly d connected components, its 
open sectors. The remedy is the following. Let us suppose that all the realizers of the 
tropical determinant of M share the same parity. In [20] this parity is called the tropical 
sign of M, denoted as tsgn M. If, however, M has realizers of both signs then tsgn M is 
set to zero. The tropical sign is essentially the same as the sign of the determinant in the 
"symmetrized min-plus algebra" [3]. The computation of the tropical sign of a matrix is 
equivalent to deciding if a directed graph has a directed cycle of even length [6l §3]. This 
latter problem is solvable in 0{d^) time [311 [27] . 

For fixed V2,V3, . . . ,Vd G T'^"^ the tropical sign gives rise to a map 

r : T'^"^ {0, ±1} : tsgn(x, V2, V3, . . . , v^) , 

where (x, 1^2, ^3, . . . , Vd) is the d x d-matrix formed of the given vectors as its rows. The 
following is a signed version of Proposition [21 

Theorem 3 ([20l Theorem 4.7 and Corollary 4.8]). Either r is constantly zero, or for e = 
±1 the preimage r~^({0, e}) is a closed tropical halfspace. Each closed tropical halfspace 
arises in this way. 

Here a closed tropical halfspace in 1"^-! is the union of at least one and at most d — 1 
closed sectors of a fixed tropical hyperplane. A closed sector is the topological closure of 
an open sector. For x G T'^"^ let 

x° := X -mm{xi,X2, . . . ,Xd) ■ Id 

be the unique representative in the coset x+Ml^ with non-negative coefficients and at least 
one zero. Then the open sectors of the tropical hyperplane H{0) are the sets Si, S2, . . . , Sd 
where 

Si := {x G T'^^^ I = and x° > for j ^ i} . 

For the closed sectors we thus have Si = {x E T'^~'^ \ x° = 0}. Clearly, —a + Si, —a + 
S2, . . . , —a + Sd are the closed sector of the tropical hyperplane H{a) and similarly for 
the closed sectors. It was mentioned in the introduction that T"^"^ is homeomorphic to R. 
One specific homeomorphism is given by the map 

{xi,X2, ■ ■ ■ ,Xd) + M.ld ^ {x2-Xi,X3-Xi,...,Xd-Xi) . 

In this paper we will often identify f*~^ with M.'^~^ via this particular map. For instance, 
its inverse translates our pictures below to T^. Moreover, it allows us to discuss matters 
of ordinary convexity in T'^"^; see [2T]. Notice that the closed and open sectors of any 
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halfspace are both tropically and ordinarily convex. For references to ordinary convexity 
and, particular, to the theory of convex polytopes see [TBI HP] , In order to avoid confusion 
which type of convexity is relevant in a particular statement we will explicity say either 
"tropical" or "ordinary" throughout. 

3. Vertices, Pseudo-vertices, and Types 

The combinatorics sneaks into the picture through the choice of a system of generators 
of a tropical polytope. Let V = {vi,V2, ■ ■ ■ ,Vn) be a finite ordered sequence of points 
in T"^"^. This induces a cell decomposition of T'^~'^ into types as follows 

typev(x) := (Ti, T2, . . . , T^) , 

where = {z G [n]\vi ^ x + Sk}- The individual sets are called type entries. A type 
vector (Ti, T2, . . . , T^) with respect to V satisfies the condition Ti U T2 U ■ ■ ■ U = [n] 
(but the converse does not hold). Sometimes it will be convenient to identify the point 
sequence V in M.'^ with the matrix (vij) whose i-th row is the point Vi. 

Proposition 4 ([121 Lemma 10]). The points of a fixed type T with respect to the generat- 
ing system V form an ordinary polyhedron Xt which is tropically convex. More precisely, 

Xt = {x G T"^^^ I Xk — Xj < Vik — Vij for all j, k G [d] and i G Tj } . 

A set in f^~^ is tropically convex if it coincides with its own tropical convex hull. The 
ordinary polyhedron Xt is bounded if and only if all entries in the type T are non-empty. 
The bounded ordinary polyhedra Xt are precisely the polytropes studied in [2T] . 

Theorem 5 ([121 Theorem 15]). The collection of ordinary polyhedra Xt, where T ranges 
over all types, is a polyhedral decomposition of T'^-'^ . The tropical polytope tconv V is 
precisely the union of the bounded types. 

Let IX be the {n + d — l)-dimensional real vector space ]R""'"'^/]R(1„, — 1^); its elements 
are the equivalence classes of pairs {y, z) + M(ln, —Id) where y G M" and z G W^. For an 
n X (i-matrix V = (vij) the polyhedron 

(2) Tv ■= {{y,z) eU\yi + Zj < Vij for all i G [n] and j G [d]} 

is unbounded. We call the polyhedron the envelope of tconv V with respect to V. By 
[T2I Lemma 22] the map IX — T*^"^ : {y, z) ^ z sends the bounded faces of to the 
bounded types induced by the generating set V . In particular, tconv V is the image of the 
union of all bounded faces of its envelope. Similarly the unbounded types correspond to 
the unbounded faces of Ty. Of particular importance to us are the vertices of 7v which 
we call the pseudo-vertices of tconv V with respect to the generating system V . What we 
just discussed is explicitly stated as Algorithm |X| below. 



input : V C T"'^^ finite 

output: pseudo-vertices of tconv V 

compute the vertices W of the envelope ^S^v via an arbitrary (dual) ordinary convex 
hull algorithm 

return image of W under projection (?/, z) ^ z 

Algorithm A: Computing the pseudo-vertices via ordinary convex hull. 

The ordinary convex hull problem asks to compute the facets of the ordinary convex hull 
of a given finite set of points in Euclidean space. Via cone polarity this is equivalent to the 
dual problem of enumerating the ordinary vertices from an ordinary halfspace description. 
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The latter also extends without changes to ordinary unbounded polyhedra which do not 
contain any affine subspace, and this is why a dual convex hull algorithm can be applied 
to as given in ([2]) in Algorithm |X1 For an overview of ordinary convex hull algorithms 
both from the theoretical and the practical point of view see [2], [19] . The complexity of 
the ordinary convex hull problem in variable dimension is not entirely settled. 

Remark 6. It is obvious to compute the types of the pseudo- vertices with respect to a 
given set of generators by checking the definition. More interesting is that it also works 
the other way around. Once the type of a pseudo-vertex is given one can determine 
its coordinates from the coordinates of the generators. To see this observe that if w 
is a pseudo- vertex with respect to V then after fixing an arbitrary coordinate, say Wi, 
the inequalities in Proposition H] degenerate to a system of equations which determine 

W2 - Wi, W3-Wi,...,Wd- Wi. 

For tconvV, just considered as a subset of T"^"^, there is a unique minimal system of 
generators with respect to inclusion; see P, Theorem 15.6] and [T^ Proposition 21]. These 
are the tropical vertices of tconv V. Among the generators a vertex w is recognized by 
the property that at least one of its type entries contains only the index of w itself; see 
Algorithm [HI This process takes 0{nd?) time. 

input : V C T'^"^ finite 

output: set of tropical vertices of tconv 

W 

foreach w & V do 

(Ti, T2, . . . , Td) ^ typey\|^|(tf;) 

if Tfc = for some k then 
L add w ioW 

return W 

Algorithm B: Computing the tropical vertices. 
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Figure 1. Tropical convex hull of four points in T^. 

Example 7. Consider the point sequence V with 

vi = (0, 3, 6),V2 = (0, 5, 2),V3 = (0, 0, 1),V4 = (1, 5, 0) 
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111 T^. This is the same example as the one considered in [5l Figure 1]. It is shown in 
Figure [TJ The tropical polygon tconv V has precisely ten pseudo- vertices. They have 
their coordinates and types as listed in Table [1] The first four pseudo-vertices are the 
given generators, and these are also the tropical vertices. The meaning of the last column 
will be explained in Section H] below. Altogether the type decomposition of tconv V with 
respect to V has ten vertices, 12 edges, and three two-dimensional faces. This information 
is collected in the so-called f -vector (10, 12,3). 

Table 1 . Pseudo- vert ices of the tropical polygon from Example [7] and Figure [H 



label 


w 


typey(w) 


facet sectors 




(0,3,6) 


1, 1, 1234 


3 


V2 


(0,5,2) 


2, 123, 24 




V3 


(0,0,1) 


123, 3, 34 




V4 


(1,5,0) 


24, 134, 4 




Wi 


(1,1,0) 


1234, 3, 4 


1 


W2 


(0,5,0) 


2, 1234, 4 


2 


W3 


(0,3,2) 


12, 13, 24 


23 


W4 


(0,3,4) 


1, 13, 234 


13 


Ws 


(1,4,0) 


124, 13, 4 




We 


(0,1,2) 


12, 3, 234 





Let Afc := conv{ei, 62, ... , 6^+1} be the (ordinary) regular fc-dimensional simplex. Here 
and in the sequel the standard basis vectors of M"' are denoted as ei, 62, . . . , e^. The 
product of two simplices A„_i x A^-i is an {n + d — 2)-dimensional ordinary polytope 
with n + d facets and nd vertices. We can read our matrix as a way of assigning the 
height Vij to the vertex (cj, ej) of A„_i x A^-i. Projecting back the lower convex hull of the 
lifted polytope yields a regular subdivision of A„_i x Ad-i, see [10] for a comprehensive 
treatment of the subject. The following is the main structural result about tropical 
polytopes. 

Theorem 8 ([121 Theorem 1]). The tropical polytope tconv with the induced structure 
as a polytopal complex by its bounded types is dual to the regular subdivision 0/ A„„i x A^„i 
induced by V . 

In view of this result the most natural version of the tropical convex hull problem 
perhaps is the one that asks to compute the pseudo- vertices together with all the bounded 
types. This can be achieved via applying a suitable modification of an algorithm of Kaibel 
and Pfetsch [23] for computing the face lattice of an ordinary polytope from its vertex 
facet incidences to the output of Algorithm IM 

For ordinary convex hull computations it is known that genericity assumptions allow for 
additional and sometimes faster algorithms. This is also the case in the tropical situation. 
The point set V C T"'""^ is sufficiently generic if no x /c-submatrix of V , viewed as 
a matrix, is tropically singular. This condition is equivalent to the property that the 
envelope is a simple polyhedron. This is the case if and only if the induced subdivision 
of the product of simplices is a triangulation. A /c-dimensional polyhedron is simple if 
each vertex is contained in exactly k facets; for details about polytopes and polyhedra 
see [10]. The maximal cells of the polytopal complex dual to the type decomposition of 
tconv are in bijection with the pseudo-vertices. If (Ti, T2, . . . , T^) is the type of the 
pseudo-vertex w then the corresponding maximal cell contains precisely those vertices 
(cj, Cj) of A„_i X Arf„i which satisfy i G Tj. 
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Example 9. In fact, the point set V in Example [7] is sufficiently generic, and hence the 
subdivision of A3 x A2 induced by is a triangulation. For instance, the tropical vertex 
Vi has the type (1, 1, 1234), and it corresponds to the maximal cell with vertices 

{(ei, ei), (ei, 62), (ei, 63), (eg, 63), (63, 63), (64, 63)} , 

a 5-simplex. See the Table [2] for the complete list of pseudo- vert ices versus maximal cells 
in the triangulation; the vertex {ei,ej) is abbreviated as ij. 

Table 2. Triangulation of A3 x A2 dual to tconv\^ decomposed into types. 



label 


maximal cell in triangulation 


min. generator of /* 


Vl 


11 12 13 23 33 43 


a;2ia;22a;3i3^323;4ia;42 


V2 


12 21 22 23 32 43 




V3 


11 21 31 32 33 43 


3^122^132^222:232^412:42 


Vi 


12 21 32 41 42 43 


2:112^132^222:232^312:33 


Wi 


11 21 31 32 41 43 


2: 1 22; 13 2^22 2:23 2^33 2:42 


W2 


12 21 22 32 42 43 


2:112^132^232:312^332:41 


W3 


11 12 21 23 32 43 


2:132^222^312:332^412:42 


W4 


11 12 23 32 33 43 


2:132^212:222:312:412:42 


W5 


11 12 21 32 41 43 


2:132:222:232:312:332:42 


We 


11 21 23 32 33 43 


2:122^132:222:312^412^42 



Block and Yu [5] investigated the relationship of tropical convexity to commutative 
algebra. We will review some of their ideas in the following. The situation is particularly 
clear if the points considered are sufficiently generic. 

Let ]K[a;ii,a;i2, . . . ,Xnd] be the polynomial ring in nd indeterminates over the field K. 
The indeterminate Xij will be assigned the weight Vij, the j-th coordinate of the point Vi 
in our sequence V. Weights of monomials are extended additively, that is, the weight of 
x° = n ^ij''' is '^CiijVij. The weight iny(/) of polynomial / is the sum of its terms of 
maximal weight. This defines a partial term ordering on ]K[xii,xi2, . . . ,2:„d]- Let J be 
the determinantal ideal generated by all 2 x 2-minors. Then 

iny(J) := (iny(/) | / G J) 

is the initial ideal of J with respect to V . 

Proposition 10 ([3 Proposition 4]). The points in V are sufficiently generic if and only 
ifinv{J) is a square-free monomial ideal. 

A square-free monomial in K[xu,xi2, ... ,Xnd] corresponds to a subset of the inde- 
terminates Xii,Xi2, ... ,Xnd and conversely. The unique minimal set of generators of a 
square-free monomial ideal is the set of minimal non-faces of a unique finite simplicial 
complex on the vertices Xu, X12, . . . , Xnd- Starting with our ideal / := iny( J) this simpli- 
cial complex is the initial complex Ay (J). The reverse construction can be applied to any 
finite simplicial complex: Its minimal non-faces generate the Stanley-Reisner ideal of the 
complex. In particular, the initial ideal / is the Stanley-Reisner ideal of the initial com- 
plex Ay (J). The complements of the maximal faces of Ay (J) generate another squarefree 
monomial ideal, the Alexander dual I*. For details see [291 Chapter 1]. The whole point 
of this discussion is that Ay(J) coincides with the triangulation of A„„i x A^.i induced 
by ^ O Lemma 5]. We arrive at Algorithm [Cl below as an alternative to Algorithm lAl 
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input : V C T'^"^ finite, in general position 
output: types of tlie pseudo- vertices of tconv V 
J ^ (2 X 2-minors of the n x (i-matrix (xj-,)) 
/ ^ iny(J) 

/* ^ Alexander dual of / 

return /* 

Algorithm C: Computing the pseudo- vertices via Alexander duality. 

Example 11. We continue our Example [71 In this case the initial ideal / := iny(J) reads 

/ = (xiiX22, 3^113^42, 3^123^213^33; 3^123^31; 2:13X21, 3:13X22, X13X31, X13X32, X13X41, X13X42, 
X22X31, X22X33, X22X41, X23X31, X23X41, X23X42, X31X42, X33X41, X33X42) . 

Notice that the generator X12X21X33 is of degree three, whence / is not homogeneous. The 
Alexander dual /* is generated by the monomials listed in the third column of Table [H 

We would like to point out that there are typos in the description of / in [5l page 
109], but the description of /* is correct. The initial complex Av{J), or equivalently the 
triangulation dual to tconv V subdivided into its bounded types, is the pure 5-dimensional 
simplicial complex whose facets are the complements of the sets corresponding to the ten 
minimal generators of /*; the maximal cells of Ay{J) are listed in the second column of 
Table [TJ It is known that the /-vector of any two triangulations of a product of simplices 
(without new vertices) is the same; in this case it reads (12,48,92,93,48, 10). 

The following is the main algebraic result on this topic. 

Theorem 12 ([5], Theorem 1]). // the points in V are sufficiently generic then tconv \^ 
supports a minimal free resolution of the ideal I*, as a polytopal complex. 

This says that the combinatorics of the type decomposition of tconv V completely con- 
trols the relations among the generators of the ideal I*. For details on cellular resolutions 
the reader is referred to [291 Chapter 4]. 

Example 13. We explain in slightly more detail where this leads to for our running 
example. Abbreviating S := ]K[xii,xi2, . . . ,X43], the cellular resolution of Theorem [T2l is 
the complex (in the homological sense) 

(3) ^S^ ^ S^^ ^ ^ ^0 

of free 5- modules. The degrees 10, 12, and 3 of the modules correspond to the /-vector of 
the type decomposition of tconv V; see Example [71 The arrows are homomorphisms of S- 
modules (to be read off the combinatorics of tconv V) with the property that concatenating 
two consecutive arrows gives the zero map; the precise maps for this specific example are 
given in [5l Example 10]. Knowing the complex ([3|) (with its maps) allows to reconstruct 
the combinatorics of the types induced by V. 

Even if one is not interested in algebraic applications. Theorem [121 opens up an addi- 
tional line of attack for the algorithmic problem to compute all the types together with 
the pseudo-vertices. This is due to the fact that a simple polyhedron, such as the en- 
velope in the sufficiently generic case, supports a unique minimal free resolution on 
the polytopal subcomplex of its bounded faces [3 Remark 7]. That is to say, in order 
to compute the types it suffices to compute any free resolution and to simplify it to a 
minimal one afterwards [5l Algorithm 2]. 

Remark 14. In ordinary convexity there are at least two algorithms to compute convex 
hulls essentially via sidedness queries: the beneath-and-beyond method (which iteratively 
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produces a triangulation) and gift wrapping [21 [19]. In general, it is not clear to what 
extent these also exist in tropical versions. However, in the planar case, that is, d = 3, the 
cyclic ordering of the tropical vertices can be computed by an output-sensitive algorithm 
in 0{nlogh) time [20l Theorem 5.3], where n is the number of generators and h is the 
number of tropical vertices, by a suitable translation from the situation in the ordinary 
case [7]. 

4. Minimal Tropical Halfspaces 

A key theorem about ordinary polytopes states that an ordinary convex polytope, that 
is, the ordinary convex hull of finitely many points in R'^ is the same as the intersection 
of finitely many affine halfspaces [IHl Theorem 2.15]. For tropical polytopes there is a 
similar statement. 

Theorem 15 ([201 Theorem 3.6]). The tropical polytopes are precisely the bounded inter- 
sections of tropical halfspaces. 

For a e T'^^i and H C [d] we use the notation (a, H) for the closed tropical halfspace 
a + Ui6-ff tropical halfspace is minimal with respect to a tropical polytope tconv V 
if it is minimal with respect to inclusion among all tropical halfspaces containing tconv V. 
This raises the question how these can be computed. We begin with a straightforward 
procedure to check if one tropical halfspace is contained in another. 

Lemma 16. Let (a, H) and {b, K) be tropical halfspaces. Then (a, H) C (6, K) if and 
only if H C K and a — b & Sk. 

Proof. The property (a, H) C (6, K) is equivalent to the following: for all x G (0, H) 
we have a + x — b G {0,K). So assume the latter. Then, as x G {0,H), this implies 
a — b & (0, K). Now we can assume that a = b without loss of generality and H (1 K is 
immediate. The converse follows in a similar way. □ 

It is a feature of tropical convexity that d of the minimal tropical halfspaces of a tropical 
polytope can be read off the generator matrix V right away. To this end we define the 
k-th corner of tconv V as 

CkiV) := i-Vi^k) Q Vi ® i-V2,k) Q V2 ® ■ ■ ■ ® {-Vn,k) Q . 

By construction it is clear that each corner is a point in tconv V. 

Lemma 17. The closed tropical halfspace CkiV) + Sk is minimal with respect to tconv V . 

Proof. By symmetry we can assume that k = 1. Then the vectors (— © t'l, (—^^2,1) © 
V2, ■ ■ . , {—Vn,i) © Vn are the Euclidean coordinate vectors of the points Vi,V2, ■ ■ ■ ,Vn, and 
their pointwise minimum c := Ci{V) is the "lower left" corner of their tropical convex 
hull. By construction the tropical halfspace c + Si contains tconv V. Suppose that c-\- §1 
is not minimal. Then there must be some other tropical halfspace w + Sk contained in 
c + Si which still contains tconv By Lemma [16] we have K = {1}. Moreover, since 
c G tconv V, we have c — w G 5*1 and so c + iSi C w + S*!, again by Lemma [T6l We conclude 
that w = c, and this completes our proof. □ 

It is now a consequence of the following proposition that the corners are pseudo- vertices. 

Proposition 18 ( [201 Proposition 3.3]). The apex of a minimal tropical halfspace with 
respect to tconv is a pseudo-vertex with respect to tconv V^. 
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The tropical halfspace Ck{V) + Sk is called the k-th cornered tropical halfspace with 
respect to tconv V. The cornered hull of tconv V is the intersection of all d cornered 
tropical halfspaces. Each tropical halfspace is tropically convex. But a tropical halfspace 
is convex in the ordinary sense if and only if it consists of a single sector. This implies 
that the cornered hull of a tropical polytope is tropically convex and also convex in the 
ordinary sense. Since it is also bounded it is a polytrope. In fact, the cornered hull of V 
is the smallest polytrope containing V. 

Example 19. For the matrix V introduced in Example [7] the three corners are the 
pseudo- vertices 

ci(\/) = ti;i = (1,1,0), C2{V) = We = (0,5,0), and csiV) = v, = (0,3,6) . 
This notation fits with Table [TJ The cornered hull is shaded lightly in Figure [H 
Lemma 20. Let {a,H) be any tropical halfspace containing V, and let (Ti,T2, . . . ,Td) = 
typey(a). Then IJjG//^i = N- 

Proof. The type entry Tj contains (the indices of) the generators in V which are contained 
in the closed sector a + Sj. Since all the generators are contained in each facet defining 
tropical halfspace we have IJje//-^i ~ [^]- '-' 

From this lemma it is clear that for a minimal tropical halfspace containing tconv V the 
type entries of a with respect to V form a set covering of [n] which is minimal with respect 
to inclusion. It may happen that one pseudo-vertex occurs as the apex of two distinct 
minimal tropical halfspaces. A tropical halfspace is locally minimal if it is minimal with 
respect to inclusion among all tropical halfspaces with a fixed apex and containing V. 

input : X eT'^ point in the boundary of a tropical polytope tconv V 

output: set of tropical halfspaces with apex x locally minimal with respect to P 

compute T = (To, . . . ,Td) = typey (x) 

compute all (with respect to inclusion) minimal set coverings of [n] hj Tq, . . . ,Ta 

^ pairs of x with all these minimal set covers 
return !K 

Algorithm D: Computing the minimal tropical halfspaces (locally). 

The input to Algorithm [E] below can be the set of pseudo- vertices with respect to any 
generating system, in particular, with respect to the tropical vertices. 

Notice that, as stated. Algorithm [D] contains the NP-complete problem of finding a 
set covering of [n] of minimal cardinality as a subproblem. The following example raises 
the question if, in order to compute all (globally) minimal tropical halfspaces it can be 
avoided to compute all locally minimal tropical halfspaces for all pseudo- vertices. This is 
not clear to the author. 

Example 21. Consider the tropical convex hull tconv(f i, ^2, fs, W4) shown in Figure [H 
The points W4 = (0,3,2) and W5 = (0,3,4) both are pseudo-vertices with respect to the 
given system of generators, see Table [H At W5 there are precisely two locally minimal 
tropical halfspaces, namely (w^, {1, 3}) and (w^, {2, 3}). The only locally minimal tropical 
halfspace at is (^4, {2, 3}). As (^4, {2, 3}) is contained in (^5, {2, 3}), the latter cannot 
be (globally) minimal. 

Both (w4, {2, 3}) and (w^, {1,3}) are (globally) minimal. Altogether there are five min- 
imal tropical halfspaces: the three remaining ones are (vi,{3}), (wi,{l}), and wq,{2}). 
For each pseudo- vertex which is the apex of a minimal tropical halfspace the corresponding 
facet sectors are listed in the final column of Table [H 
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input : W G T'^ pseudo- vertices of tconv V 

output: set of minimal tropical halfspaces with respect to tconv 

foreach x G W do 

J-C' ^ locally minimal tropical halfspaces at x according to Algorithm [D] 
foreach (x, H') e "K' do 
foreach (y, iJ) G IK do 

if (x, H') is contained in (y, H) then 
I remove (?/, H) from % 
else 

if (y, H) is contained in (x, H') then 
|_ remove (x, H') from "K' 

L !H < — IK U CK' 

return J{ 

Algorithm E: Computing the minimal tropical halfspaces (globally). 

5. Matroid Subdivisions and Tree Arrangements 

It turns out that tropical polytopes or, dually, regular subdivisions of products of ordi- 
nary simplices carry information which can be, in a way, extended to matroid subdivisions 
of hypersimplices. These are known to be relevant for studying questions on the Grass- 
mannians of rf-planes in n-space over some field IK [2U [371 [U] • 

An ordinary polytope in whose vertices are 0/1-vectors is a matroid polytope if each 
edge is parallel to a difference of basis vectors Cj — ej for some distinct i and j in [m] . A 
matroid of rank r on the set [m] is a set of r-element subsets of [m] whose characteristic 
vectors are the vertices of a matroid polytope. The elements of a matroid are its bases. 
It is a result of Gel'fand, Goresky, MacPherson, and Serganova [15] that this definition 
describes the same kind of objects as more standard ones [39]. Each m x n- matrix M 
(with coefficients in an arbitrary field K) gives rise to a matroid M of rank r on the set 
[m] as follows: the columns of M are indexed by [m], and r is the rank of M. Those 
0/1-vectors which correspond to bases of the column space of M are the bases of M. We 
denote the matroid polytope of a matroid M as P^vt- The {r,m)-hypersimplex A(r, m) is 
the matroid polytope of the uniform matroid of rank r on the set [m] , that is, the matroid 
whose set of bases is (''^'), by which we denote the set of all r-element subsets of [m]. 

The relationship to tropical convexity works as follows. By the Separation Theorem 
from ordinary convexity [THl §2.2] each vertex v of an ordinary polytope P can be separated 
from the other vertices by an affine hyperplane H. The intersection of P with H is again 
an ordinary polytope, and its combinatorial type does not depend on H; this is the vertex 
figure P/v of P with respect to v. Moreover, whenever we have a lifting function A on 
P inducing a regular subdivision, then A induces a lifting function of each vertex figure 
P/v, and this way we obtain a regular subdivision of P/v. In the case of the {d,n + d)- 
hypersimplex each vertex figure is isomorphic to the product of simplices Ad-i x A.„_i. 
Hence regular subdivisions of A{d,n + d) yield configurations of n points in 1"^"^ at each 
vertex of A{d,n + d). The situation for the hypersimplices is special in that the converse 
also holds, that is, each regular subdivision of a vertex figure can be lifted to a (particularly 
interesting kind of) regular subdivision of A{d,n + d). 

Proposition 22 ([231 Corollary 1.4.14]). Each configuration of n points in T'^"! can be 
lifted to a regular matroid decomposition of A{d,n + d). 



12 



MICHAEL JOSWIG 



A matroid subdivision of a matroid polytope is a polytopal subdivision with the property 
that each cell is a matroid polytope. Kapranov's original proof [21] uses non-trivial 
methods from algebraic geometry. Here we give an elementary proof, which makes use of 
the techniques developed in [17]. 

Proof. Let V be the d x n-matrix which has the n given points as its columns, and let V 
be the d x n-matrix produced by concatenating V with the tropical identity matrix of size 
d X d. The tropical identity matrix has zeros on the diagonal and infinity off the diagonal. 
Each d-subset S oi[n + d\ defines d columns of V and this way a x rf-submatrix Vs- It 
is easy to see that the map 

Ti : S ^ tdet Vs 

is a finite tropical Pliicker vector. A finite tropical Pliicker vector is a map from the set 
('"d*^') reals satisfying the following property: For each subset T of [n + d] with 

d — 2 elements the minimum 

min{7r(Tzj) + 7r(TA;/), 7T{Tik) + 7r(Tj7), 7r(Tz/) + 7r{Tjk)} 

is attained at least twice, where i,j, k, I are the pairwise distinct elements of [n + d]\T 
and Tij is short for TU{i, j}. Clearly, the map vr is a lifting function on the hypersimplex 
A{d, n + d), and thus it induces a regular subdivision. That the tropical Pliicker vectors, 
in fact, induce matroid decompositions is a known fact [Ml §1.2], [351 Proposition 2.2]. 

One can verify directly that the point configuration V in 1"^-! is (in the sense of 
Theorem [8]) dual to the regular subdivision induced by vr at the vertex figure of A{d,n + d) 
at the vertex e„+i + e„+2 + ■ ■ ■ + Cn+d- CH 

Notice that our computation in the proof above makes use of "cxd" as a coordinate, that 
is, from now on we work in the tropical projective space 

TP'^^^ := ((RU{oo})"'\{(oo,oo,...,oo)})/Mlrf. 

The tropical projective space TP'^^^ is a compactification of the tropical torus T^~^ with 
boundary. In fact, T'^"^ should be seen as an open ordinary regular simplex of infinite 
size which is compactified by TP'^"^ in the natural way. 

Example 23. Starting out with the point configuration V from Example [7] with d = 3 
and n = 4 we obtain 

/O 1 cx) oo\ 
(4) y= 3505cx)0oo . 

\6 2 1 oo CX) y 

The resulting tropical Pliicker vector is given in Tabled and the induced matroid subdivi- 
sion of A{d,n + d) is listed in Table HI For each of the ten maximal cells the corresponding 
set of bases is listed. These maximal cells are in bijection with the pseudo-vertices of 
tconv V with respect to V as given in Tables [T] and [2J 

Table 3. Tropical Pliicker vector vr of from Equation (jl]). 



123:2 124:3 125:5 126:2 127:3 134:0 135:4 136:1 137:0 145:3 146:0 147:4 
156:6 157:3 167:0 234:0 235:2 236:1 237:0 245:5 246:0 247:5 256:2 257:5 
267:0 345:0 346:0 347:1 356:1 357:0 367:0 456:0 457:5 467:1 567:0 



The hypersimplex A{d, n+d) is the set of points {xi, X2, ■ ■ ■ , Xn+d) £ M""'"'^ which satisfies 
the 2{n + d) linear inequalities < Xi < 1 and the linear equation = c?. The given 
inequalities are all facet defining. The facet defined by Xj = is the i-th deletion facet, 
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Table 4. Matroid subdivision of A(3, 7) induced by the tropical Pliicker 
vector vr in Table [31 



label 








matroid bases 




Izo 


IZD 


135 


'\ ofi I/If; ic^/^ ic^T i/^v occ/^ oc^^ r^p.'~7 

lab 145 14d iob ibi Ibi zoo oob 4od ob/ 


V2 


lz4 


IzO 


Iz ( 


1 /I 1 f:'7 oQ/1 OQf^ OQV O/If; 0/l^i of;^^ oc^v o^^'7 Q/If; qf;'? /Ic;/^ c^c'v 

140 Lo { zo4 Zoo zo/ z4o z4b zob zo^ zb/ o4o 60 ( 4ob ob/ 




1 Q /I 

lo4 


lob 


13 1 


i4b lb I 264: zob zo/ z4b zb/ o4o o4b oob 00/ 6b ( 4ob obi 




lz4 


Iz / 


145 


'\ A'7 1 F;'7 OO/I 007 O/lf^ O/IV QAV:. QAT AKfi A^'7 AP.'7 ^P.'y 

14/ ibi zo4 zo/ z4b Z4( Ibi o4o o4/^ 00/ 4ob 4b i 4b/ ob/ 


Wi 


134 


137 


146 


167 234 237 246 267 345 346 347 357 367 456 467 567 


W2 


124 


127 


145 


157 234 237 245 246 247 257 267 345 357 456 457 567 


W3 


123 


124 


125 


126 127 134 137 145 146 157 167 234 235 237 246 256 267 345 357 456 567 


1^4 


123 


125 


126 


134 135 136 137 145 146 157 167 235 256 345 356 357 456 567 


W5 


124 


127 


134 


137 145 146 147 157 167 234 237 246 267 345 347 357 456 467 567 


We 


123 


126 


134 


136 137 146 167 234 235 236 237 246 256 267 345 356 357 456 567 



and the facet Xi = 1 is the i-th contraction facet of A((i, n + d). Each deletion facet is 
isomorphic to A((i, n+rf— 1), and each contraction facet is isomorphic to A((i— 1, n+rf— 1). 
Recursively, all faces of hypersimplices are hypersimplices. Now Pjvc is a subpolytope of 
A((i, n + d) whenever M is a matroid of rank d on the set [n + d], and the intersection 
of Pjvt with a facet is again a matroid polytope, a deletion or a contraction depending on 
the type of the facet. This implies that each matroid subdivision of A{d,n + d) induces 
a matroid subdivision on each facet. 

Table 5. Restriction vti of the tropical Pliicker vector vr from Table [3] to 
the first contraction facet. 



23:2 24:3 25:5 26:2 27:3 34:0 35:4 36:1 37:0 45:3 46:0 47:4 56:6 57:3 67:0 



We want to look at the matroid subdivisions of A{d,n + d) for low rank d. The 
hypersimplex A(l,n + 1) is an n-dimensional simplex without any non-trivial (matroid) 
subdivisions. For d = 2 the situation is more interesting: Each matroid subdivision of 
A(2,n + 2) is dual to a tree with n leaves [24l §1.3]. For c? = 3 we have the following 
result. 

Theorem 24 ([TTl Theorem 4.4]). The regular matroid subdivisions ci/A(3,n + 3) bi- 
jectively correspond to the equivalence classes of arrangements of n + 3 metric trees with 
n + 2 labeled leaves. 

It is easy to see where the n + 3 trees come from: They are dual to the matroid subdi- 
visions induced on the n + 3 contraction facets of A(3, n + 3). A sequence {Si, S2, ■ ■ ■ , Sn) 
is an arrangement of metric trees if 6i is a tree metric on the set [n] \ {i} satisfying 

Si{j,k) = 6j{k,i) = 6kii,j) 

for any three distinct i,j,k G [n]. Arrangements of metric trees are equivalent if they 
induce the same arrangement of abstract trees, which roughly means that their intersection 
patterns are the same; see [17] for details. 

Example 25. Restricting the tropical Pliicker vector vr on A(3, 7) from Example [231 to 
the first contraction facet gives a tropical Pliicker vector vri on the second hypersimplex 
A(2, 6) (with coordinate directions labeled 2, 3, . . . , 7) shown in Table[5l Since tti happens 
to take values between and 6 we can set 6i{S) := 3- (7ri(^)/6) for S G (^^'^^■■'^^) , which 
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yields a metric on the set {2, 3, . 
have 



Si 



7}. Written as one half of a symmetric matrix we 



( 


8/3 5/2 


13/0 


8/3 


5/2\ 




U o 


7/3 


1 / /O 


Q 
O 







5/2 


3 


7/3 









2 



5/2 
3 


\ 








/ 



Notice that the matroid subdivision of the first contraction facet induced by vr does not 
change if we add multiples of (1, 1, ... , 1) to vr or scale it by any positive real number. 
The rescaling chosen above makes sure that Si takes values between two and three only. 
This way the triangle inequality becomes valid automatically, which is why Si is indeed a 
metric. Moreover, the metric Si is tree-like, and the corresponding metric tree is shown 
in Figure [3 For instance, one can check that 5i(2,3) = 8/3 = 13/12 + 1/6 + 17/12. 
The edge lengths on the tree can be computed via the Split Decomposition Theorem of 
Bandelt and Dress [1]. 



3/4 



1/6 



1/6 



1/4 



7/6 



17/12 



13/12 



Figure 2. Caterpillar tree with six labeled leaves and edge lengths corre- 
sponding to the metric Si. 

So far we looked at a single tree only. Now we want to explore how the n + 3 trees, 
one for each contraction facet, interact. In the plane a tropical line is the same as a 
hyperplane. That is to say, the tropical hyperplane/line 

H{a) = -a+ (R>oei U M>oe2 U M>oe3) 

is the union of the infinite rays into the three coordinate directions emanating from the 
apex —a. In particular, it can be seen as a metric tree with one trivalent internal node 
and three edges of infinite length. The taxa sit at the endpoints of those rays in the 
compactification TP^. The mirror image of H{a) at its apex is its dual line 

H*{a) = -a - (M>oei U M>oe2 U M>oe3) . 

Notice that H*{a) is a tropical line if one takes "max" as the tropical addition rather than 
"min" . Any two tropical lines either share an infinite portion of some ray or they meet in 
a single point, and the same holds for their duals. The arrangement of dual tropical lines 
induced by a point sequence ^ in T^, denoted as H*{—V), is the sequence of tropical lines 



with apices fi, f2. 



In order to make the connection to matroid decompositions of 



hypersimplices we can instead also look at the compactified version H*(—V) which is the 
sequence of lines in TP^ with the n-l-3 apices Vi,V2, ■ ■ ■ , Vn, (0, oo, oo), (oo, 0, oo), (oo, oo, 0). 
Now each dual tropical line H*{—Vi), for 1 < i < n + 3, gives rise to a labeled tree by 
recording the intersection pattern with the other trees. 

Example 26. The tree arrangement H*{—V) and its compactification H*{—V) arising 
from our running example matrix V introduced in Example [7] are shown in Figure [3]to the 
left. The tree with apex = (0, oo, oo) occurs as the dashed line on the top and on the 
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Figure 3. Arrangement of seven trees, three of which cover the boundary 
of TP^ 

right of the rectangular section of T^, the tree corresponding to Vq = (oo, 0, oo) is to the 
left, and the tree corresponding to vj = (oo, oo, 0) is at the bottom. For instance, the tree 
corresponding to Vi receives the labeling shown in Figure [2l The subsequent Theorem [27] 
can be visually verified for our example by comparing Figures [1] and [31 

Theorem 27 ([H Theorem 6.3]). Suppose that V is a sufficiently generic sequence of 
points in T^. Then the polyhedral subdivision of induced by the tree arrangement 
H*{—V) coincides with the type decomposition induced by V . 

The abstract arrangement of trees induced by the matroid decomposition of A(3, n + 
3) which in turn is induced by the point sequence V is precisely the compactified tree 
arrangement II* {—V). 

Remark 28. The Cayley trick from polyhedral combinatorics allows to view triangula- 
tions of A„„i X Arf_i as lozenge tilings of the dilated simplex nAd-i, see [32]. For our 
Example [26] the corresponding lozenge tiling of 4A2 is shown in Figure [3] to the right. The 
tree arrangement II*{—V) partitions the dual graph of the tiling. 

6. Computational Experiments with polymake 

Most algorithms mentioned above are implemented in polymake, version 2.9.4, which 
can be downloaded from www . polymake . de . After starting the program on the command 
line (with the application tropical as default) an interactive shell starts which receives 
commands in a Perl dialect. For instance, to visualize our running Example [7] it suffices 
to say: 

> $p = new TropicalPolytope<Ratioiial> (POINTS=>[ [0,3, 6] , [0,5,2] , [0,0,1] , [1,5,0]]) ; 

> $p->VISUAL; 

We want to show how polymake can be used for the investigation of tropical polytopes. 
For instance we can compute the type of a cell in the type decomposition as follows. This 
is not a standard function, which is why it requires a sequence of commands. Again the 
pseudo-vertices and the types always refer to the fixed sequence of generators that we 
started out with above. 

> print @{rows_labeled($p->PSEUDOVERTICES)}; 
0:0 5 
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1:0 4 -1 
2:0 3 2 
3:0 3 4 
4:0 3 -1 
5:0 1 
6:0 -1 
7:0 1 2 
8:0 5 2 
9:0 3 6 

> print $p->ENVELOPE->BOUNDED_COMPLEX; 
{01248} 

{2 3 7} 
{24567} 
{3 9} 

The first command lists the pseudo- vertices (in no particular order), and the function 
rows_labeled is responsible for making the implicit numbering explicit which arises from 
the order. The next command lists all the bounded cells of the type decomposition. Each 
row corresponds to one maximal cell, namely the (tropical or ordinary) convex hull of the 
pseudo-vertices with the given indices. We obtain two pentagons, one triangle, and one 
edge. Now we want to compute the type of the first cell (numbered '0'): 

> $indices = $p->ENVELOPE->BOUNDED_CDMPLEX-> [0] ; 

This is just the sequence of indices of the pseudo-vertices defining this cell. In the 
following we produce an ordinary polytope which has these points as its vertices. One 
minor technical complication arises from the fact that the leading '0' used for coordinate 
homogenization in the tropical world must be replaced by a '1' which makes the ordinary 
polytope homogeneous. 

> $vertices = $p->PSEUDOVERTICES->minor ($indices , range (1 , $p->AMBIENT_DIM) ) ; 

> $n_vertices = scalar (@{$vertices}) ; 

> $all_ones = new Vector<Rational> ( [ (l)x$n_vertices ]); 

> $cell = new Polytope<Rational> (VERTICES => ($all_ones | $vertices) ) ; 

> print $cell->VERTICES; 
15 

14-1 
13 2 
13-1 
15 2 

Now the type to be computed is just the type of any relatively interior point of this 
ordinary polytope called $cell. 

> $point = $cell->REL_INT_POINT; 

> print $point ; 
1 4 2/5 

This point must be translated back into the tropical world, and then we can call a 
function to compute its type, which turns out to be (2,13,4). Notice that polymake 
starts the numbering from '0', and hence each index in the output is shifted by one. 

> $cell_point = poly2trop(new Polytope<Rational>(POINTS=>[$point] ))->POINTS; 

> print types ($cell_point, $p->PDINTS) ; 
{1} 

{0 2} 
{3} 



TROPICAL CONVEX HULL COMPUTATIONS 



17 



It should be mentioned that there is also the Maxplus toolbox for Scilab for com- 
putations in tropical convexity [25] . 
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